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Abstract 

In this paper a new approach to solving the 2D and 3D Ising models in 
external magnetic field H ^ is developed. The general formalism for the 
approach to the problem is presented on the example of the 2D Ising model in 
the external magnetic field. The paper presents a new method obtaining the 
Onsager solution and computations of asymptotic forms of low-temperature 
free energy for the 2D Ising model in the external magnetic field (H). The 
free energy in the limiting case of the magnetic field tending to zero (H — ► 
0, N,M — > oo) at arbitrary temperature is also considered (T ^ 0). 
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I. INTRODUCTION 



We will briefly describe the well known model of a magnetic containing(?) variety of 
spins situated on the vertices of a crystalline lattice. The spin ctt k can be "up" (cr*. = 1) or 
"down" (o"fc = — 1). A microscopic state of the system is characterized by orientations of all 
the spins. Energy E{a} of the microscopic state {a} is composed of two contributions, one 
from the exchange interactions of the spins and described by the interaction constant J k i, 
and the second from the interaction of the spins with the external magnetic field (if): 

E{a} = -J2 J ki a k°i ~ H^a k , (1.1) 

ki k 

where summation is taken over all sites of the lattice. The key problem is calculation of the 
statistical sum: 



Z=J2*M-PE{a}), P = t^i (1-2) 



where T denotes temperature and kg - the Boltzman constant. 

The model described above was introduced by W.Lenz in 1920 |TJ, and for the one- 
dimensional case was investigated by E.Ising in 1925 0. The first exact solution of the 
statistical mechanical problem for the 2D, (H = 0) case was found by L.Onsager in 1944 
0. We use the standard name, the Ising model. 

The solution given by Onsager strongly influenced the development of all of statistical 
physics, and in particular of the theory of phase transitions. It was shown for the first time 
that exact calculation of the free energy leads to an evidence that thermodynamic quantities 
behave in the vicinity of the phase transition in a way which is essentially different from 
that in the approximate models, like e.g. the mean-field theory. The result for spontaneous 
magnetization 97to in the model was presented by Onsager at the conference in Florence in 
1949 |4j], i.e. 5 years after the successful derivation of the expression for free energy. The 
first published derivation for ffl was given by Yang |5j , are recently, alternative derivations 
have been published both for the free energy and 9Jto • 

In spite of its simplicity, the Ising model is not only non-trivial in higher dimensions 
(d > 2), but also it has rich structure. By this we mean not only its connection with other 
models (for example with the lattice gas models, binary alloys, some models in quantum 
field theory [|7|, |TU| etc.), and wide application in numerous domains of statistical physics, 
but also its role as a generator of new ideas and tools, which find its use in various areas 
of physics and mathematics. There are sufficiently many examples of such applications and 



we will not discuss them here (some examples can be found in the monograph 11], where 
stochastic Ising models are considered, and also their connection with Markov processes with 
local interactions). We would like to stress that this rich structure of the Ising model has 
maintained a high level of interest in this problem among physicists and mathematicians. 

In this paper we present a new approach to the Ising problem in external magnetic field 
(if), with the nearest-neighbour interaction on the square lattice. In connection with that 
we would like to mention the paper by Schultz, Mattis and Lieb [Bj, who applied it to 
solve the 2D Ising model without an external magnetic field. To calculate 9Jto they used 
a method based on a transfer-matrix using a transformation to a fermionic representation. 
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This deep, clear and logically closed paper influenced strongly the author and moved him to 
search for the solution of the problem in external magnetic field. The fundamental idea of 
the approach of the authors of the paper |J is transition to a fermionic representation (the 
transfer-matrix method was essentially used already in the paper by Onsager ||), and this 
can be treated in a sense as a problem of interacting fermions on the one-dimensional lattice. 
In this paper we use essentially the same idea. The difference is the fermionic representation 
is introduced not on a ID lattice (where the T-matrix is expressed in terms of the Fermi 
creation and annihilation operators (c£, c n ), ||) but on a two-dimensional lattice with the 
doubly indexed Fermi creation - annihilation operators (c^ m ,c nm ), JTj]. 



II. FORMULATION OF THE PROBLEM 



Let us consider the square lattice composed of M columns and N rows, on the vertices 
of which the quantities a nm taking one of the two values ±1 are defined. We will call the 
quantities the Ising " spins" . The multiindex nm numbers the sites of the lattice, where 
n numbers a row, and m numbers a column. The Ising model with the nearest-neighbour 
interaction in external magnetic field is given by the Hamiltonian of the form: 

7~C — ~ J2 E ffnm&n+l,m ~ J\ E ^nm°'n,m+l ~~ H <7 nm , (2-1) 
nm nm nm 

which takes into account anisotropy in the interaction ( J 12 > 0) between nearest neighbours, 
and also the interaction of the spins a nm with external magnetic field H, directed "up" 
(d nm = +1). The essential problem is calculation of the statistical sum for the system: 

Z(h)= £ ... £ ex P (-/3W) 

CTH=±1 <TjVM=±l 



ex p 

(fnm = ±l) 



NM 



nm&n+\,m ~l~ K\O nm (J n m +\ + h(7 nm ) 

n,m=l 



where 



#1,2 = (3Ji,2, h = 0H, (3 = l/k B T. 



(2.2) 



(2.3) 



Periodic boundary conditions are introduced for the variables u nm . Let us mention here that 
the statistical sum (2.2) is symmetric with respect to the change h — > —h where h is defined 
above (2.3). 

As is known 0, the statistical sum for the 2D Ising model in external field (H) in the 
representation of second quantization can be written in the form: 



Z = Tr(V) N = Tr(V x V 2 V h 



\N 



(2.4) 



where the operators Vi, expressed in terms of the Fermi creation and annihilation operators 



[ c mi c m), 



are of the form: 



Vi = (2sinh2A' 1 ) M/2 exp 



M 



-2^^(4^-1/2) 



m=l 



(2.5) 
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V 2 = exp < K 2 



"M-l 

£ 

_m=l 



(<",'„ - Cm)(c m+ i + C m+ i) — ( — 1) M ( 



4 _ 



'M 



a, 



(2.6) 



f M 


m— 1 






Z7T ^ 4 C P 


( C m + C m) J 


[ m=l 


P =1 





(2.7) 



where -Kj, (j = 1,2,) and h are defined above (2.3) and M = Ysl c ln c m is the operator of 
the total number of particles and K* and K\ are connected by the following formulas: 

tanh(ATi) = exp(-2fT*), or sinh 2K X sinh 2^* = 1. (2.8) 

One can see that the operator 14 in the second quantization representation, that describes 
interaction of the spins with external magnetic field, has rather complicated structure. It is 
easy to see that this operator does not commute with the operator P = (— 1) M . As a result 
the operator V 2 has also not a very tractable form, i.e. it has not the needed translational 
symmetry (2.6). More exactly, although the operators V\ and V 2 commute with the operator 
P, the operator V (2.4) does not commute with the operator P, i.e. [P, V]_ 7^ 0, because 
[P, V/j]_ 7^ 0. Therefore, we can not divide all states of the operator V = V^V 2 Vh into 
eigenstates of the operator P with eigenvalues A = ±1, and this leads to nonconservation 
of the states with even and odd numbers of fermions (for details see ||). Namely this is 
the fundamental reason which stops solving the problem under consideration within this 



formalism. Nevertheless, the author does not share Ziman's pessimism |T3[ which is based 
on some misunderstanding, because he considers actually the approach of the authors of the 
paper 0, but in the end he writes about limitations of the method of Onsager ||. In 
fact Onsager in his approach does not apply the field theoretic language of the creation and 
annihilation operators as it is in the approach of the authors of the paper 0. the method 
of Onsager ||[14[] really shows some limitations when one tries to apply it to solving the 2D 
Ising model in external magnetic field, or for solving the 3-D Ising model. On the other hand, 
completely different state of affairs we have for the approach of the authors of the paper , 
where in all its beauty the field theoretic language of the method of second quantization is 
used. The approach of the authors of the paper || allows for generalizations. We intend 
to present one of such generalizations in this and the following papers devoted to the Ising 
problem. 

Coming back to the difficulties mentioned above which are connected with the operator 
Vh, (2.7), it is now clear that to overcome the troubles within the approach 0, one should 
find an appropriate method of substituting the operator Vh (2.7) with another one which 
would be equivalent to the former in the sense of correct counting of the interaction of 
external magnetic field with the spins of the system. Namely, as it could be easily seen, the 
only contribution to Z (2.4) from the operator Vh comes, in the representation of second 
quantization, from the "even" part with respect to operators c^Cm of the operator Vh. In 
principle such transformation could be always done. However, in practice this task seems to 
be hopeless, and the direct method of calculation of the commutators used by Onsager for 
solution of the problem without external field here is simply inapplicable. We believe there 
is not an effective method to do that at least if one stays in the space of given dimension 
(d = 2 for the initial variables cx nm , and d = 1 for the variables in the representation of second 
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quantization c^ m and c m ). Nevertheless, as we will show below, there is an effective method of 
transforming the magnetic operator Vh, (2.7), after which the transformed operator allows for 
the Fourier transform of the operator V (2.4). The idea consists of formulating the problem 
in the space of higher dimension than the former one, then to pass to the representation of 
second quantization with the operator V, and afterwards to perform a limit transition with 
respect to one of the interaction constants, by going with it to zero. Having this done a 
possibility appears for effective rebuilding of the operator which is responsible for interaction 
of the spins of the system with external magnetic field. Below we will shortly present this 
approach on an example of the one-dimensional Ising model which is then applied to solving 
our basic problem. 



III. ONE-DIMENSIONAL ISING MODEL 

In the beginning of the consideration of the ID Ising model we have already the complete 
set of formulas (2.4—8). To apply them to the ID Ising model one should take simply K\ — 
and N = 1. Then, after not complicated transformations, taking into account the expressions 
(2.8), one can write the following formula for the statistical sum (2.4) [Z(K 1 = 0) = Z*\: 

Z* = Tr(V*V 2 V h ), (3.1) 

where the operators V 2 and 14 are defined above (2.6 — 7), and the operator V* is of the 
form: 



M r 



y*= n 



m=l 



1 + ("1) C 



(3.2) 



Introducing in an appropriate manner the basis in the representation of occupation num- 
bers [(finite dimensional Fock space): |0 > is the vacuum state, c m |0 >= 0; cJjO > is a 
one-particle state (m = 1,2,3, ...) etc.], and calculating the trace in (3.1) we get 

Z* = E < W?V 2 V h )\l >= 2 M < 0\(V 2 V h )\0 >, (3.3) 

all(l) 

where on the left hand side of (3.3) the summation is over all states \l >. It is easy to 
see that all the matrix elements < /|(...)|/ > in (3.3) are equal to zero thanks to the phase 
factors (— l) c ™ Cm entering the operator V*, with the exception of the vacuum matrix element 
< 1 ( — ) 1 >. For this matrix element contribution from the operator (3.2) is equal simply 
2 M . From this we obtain the right hand side of the equality (3.3). 

Let us mention now that the operators Vh, (2.7), can be represented in the form: 

M 

V h = cosh M (/i) J] [l + V m (4 + c m ) tanh h] , (3.4) 

m=l 

where the phase factor ip m is defined in an obvious way (2.7) and we applied the following 
identity 

exp(pt) = cosht + psinht, p 2 = 1. 
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Now, "dragging" the operator 14 , (3.4) through the ket- vector |0 >, after a number of 
transformations, we obtain the following representation for 14 | >: 



M 



V h \0>= cosh M (/i) ] [ e ac ™ I >, a = tanh(/i). 



(3.5) 



m=l 



Deriving the formula (3.5) we dragged all phase factors ip m entering the operator 14 through 
the vacuum state |0 > and omitted the annihilation operators c m , because c m | >= 0. We 
will omit below, for brevity, the ket- vector |0 >. This should not lead to misunderstandings. 
Further we not that the operators c+ and c\ commute with the commutator [c^, c\] = 2^^. 
As a result, using the Hausdorff- Baker formula (ct,/3 = const): 



exp (arc) exp(/3y) = expire + f3y + (afl/2)[x, y}), 



(3.6) 



after (M — 1) times application of this formula (3.6) to the operator (3.5), this operator can 
be represented in the form: 



V h = cosh M (/i) exp 



M 






exp 


m=l 





a 



M M-m 
m=l p=l 



(3.7) 



where a is defined above (3.5). Since all terms in the operator V2, (2.6) contain bilinear 
products of the Fermi operators, and the following equality is satisfied 



/ m \ M 
exp la c m) = 1 + a Yl c ln, 

\ m=l / m=l 



it is easy to see that in the pairings < | (...) | > the components linear in a m give null 
contribution. As a result, we can write the following expression (Vh— >V h *) for the operator 
14,(3.7): 



14* = cosh M (/i)exp 



M M-m 
m=l p=l 



(3.8) 



Now one can easily see tat the operator P = (— 1) M , (M = J2v" c ln c m) commutes with the 
operators Vi and 14*, and, as a consequence, the states with even or odd numbers of fermions 
are conserved. Hence, the statistical sum Z* , (.3), can be represented in the form: 



Z* = 2 M < I {VfV£) I >, 



where 



V? 



exp 



M 



C m c m)( c m +l ~\~ C. 



m+lj 



m=l 



(3.9) 



(3.10) 



and the (+) sign in V 2 + corresponds to the even states, to which are assigned the antiperiodic 
boundary conditions and the (— ) sign to the odd states, to which are assigned the periodic 
boundary conditions. 
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Passing in a standard way to the momentum representation 

exp(— in/A) 

q 



tqill 



we obtain after some simple transformations on Z*, (3.9) the expression: 



Z* = [2cosh(/i)] M < 



( n v*^) 

0<q<TV 



>. 



(3.11) 



where 



771,77* J sing 



)sing I 



V 2q = exp \2K 2 [(r/lriq + r]l q r]_ q )cos q + (r] q r]^ 

1 + cosg t t.ffs.nf \ 

-V-qVl + /(?) + /(-?) 



exp 



sm q 



in which the terms /(±g) in the expression for V£ : 



/(?) 



t t 



2 sing 

and in the case of antiperiodic boundary conditions should be omitted. 

Finally, calculating the vacuum matrix element for fixed q, after some not complicated 
transformations, we get for Z* (3.11) in the case of even states the expression (ZV): 



Zl 



[2cosh(/i)] M Jl [cosh 2K 2 -smh2K 2 cos q + a 2 sinh 21^(1 + cos q)} 

0<q<TT 

1/2 



M 



[2cosh(/i)coshir 2 ] M J] 

m=l . 

„2 



1 + ^2+ 2z 2 y - 2z 2 (l - y) cos[ 



7r(2m - 1)- 

M " 



(3.12) 



where z 2 = tanh and y = a z = tanh 2 h. Obviously, for iV noninteracting Ising models 
in external magnetic field the statistical sum W(h) is equal to the N — th power of the 
expression (3.12), i.e. W{h) = [Z+] N . In the case of odd states, as one can easily show the 
following equality is satisfied: 



Z* = 2Z* 



(3.13) 



Let us note here that the representation (3.12) unexpectedly finds an application in graph 
theory. Namely, with help of the representation (3.12) one can calculate the generating 
function for Hamilton cycles on the simple rectangular lattice (N x M), jH^ . 

Finally, we obtain the following expression for free energy per spin in the thermodynamic 
limit: 



3F = lim — In Z* = In 

M-»oo M 



e K2 cosh h + (e 2 ^ 2 sinh 2 h + e~ 2K2 ) 1/2 



(3.14) 



i.e. the known classic expression We paid so much attention to the ID Ising model 

because we wanted to show on the first place to show effectiveness of the proposed method 
of transformation of the magnetic operator V h (2.7) to its equivalent (3.8). Additionally, as 
we mentioned above, a bit different representation of the statistical sum for the ID Ising 
model (3.11) finds its application in graph theory Finally, this will help us to save 

time and place considerably when we will discuss the 2D and 3D Ising models in external 
magnetic field. 
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IV. TRANSFER-MATRIX 



In this section we will consider shortly the representation of the statistical sum for the 
3D Ising model in external magnetic field H, applying for this purpose the well known 
transfer-mart ix method |7], || [16|, 17] . An exhaustive and outstanding presentation of the 



method the reader can find in the monographs ]7|, § , where are also shown other necessary 
pieces of knowledge on application of this method. 

Let us consider a simple cubic lattice consisting of N rows, M columns, and K planes, in 
the sites of which the "spins" <J nm k are situated, which take on two values: <J nm k = ±1. The 
Hamiltonian for the 3D Ising model in external magnetic field H with the nearest-neighbour 
interaction is given in the form: 

NMK 

H = ~ E ( , A cr nmfc cr n+l,mfe + J2 (J nmk (J n,m+l,k + ■^3 cr nmfc cr nm,fe+l + Ha nmk ) , (4-1) 
(n,m,fc)=l 

where the multiindex (nmk) numbers the sites of the simple cubic lattice (N xMx K) , and H 
is the external magnetic field directed "upwards" {o n mk = +1)- The constants (Jj > 0) take 
into account anisotropy of the interaction of the Ising spins. There are periodic boundary 
conditions imposed, as usual, on the variables a nm k- The statistical sum for the system 
Z 3 (h) we write in the form: 



Uh)= E • E < 

0"111=±1 (TjvMif=±l 



-PH 



E ex p 

{ CT nmfc=±l} 



E {■Kl&nmk <7 n+l,mk + ^2 cr nmk cr n,m+l,k + ^3 cr nmfc cr nm,fc+l + ha 



nmk ) 



.nmk 



, (4-2) 



where the quantities Ki and h are defined as above (2.3), [here and everywhere below sum- 
mation over nmk (or nm) and also multiplication over nm will mean summation or multi- 
plication over the full set of integer numbers from 1 to N, M and K over each corresponding 
index, respectively]. 

In analogy to the two-dimensional case, it is convenient to introduce the notion of /c-layer 
which is understood as a set of Ising spins in all sites of a fc-layer: 

«{«} = a k = {a nmk }, k - fixed. 

Then summation in (2.4) can be conveniently executed over the layers a k , after writing the 
expression for Z 3 (h) in the form: 



K 



W = E-E«p E 



E(-^ l(T ™+ 1 . mfc + K2&n,m+l,k + K 3 a nnltk+1 + h)(T, 



nmk 



EV^ rp{(J nml} rp{(? nm2} 

ai a K 



{o"nm,if + l}' 



(4.3) 



where 



S 



Wnm,k + l} " 



^2(Ki<T n +i,mk + K 2 a n , m+ i >k + h)a, 



nrnk 



exp 



• (4-4) 



We will impose now periodic boundary conditions on the indices n, m, and k, taking 

(4.5) 

As a consequence of what was stated above and of the conditions (4.5) we can write Z 3 (h) 
in the form 

Z 3 {h) = Tr(T) K , (4.6) 

where T is the transfer-matrix, matrix elements of which are described by equalities (4.4). 
Matrix elements of the transfer-matrix of the layer-layer Ising model can be written in a bit 
different form |7[], but all these representations are in fact equivalent. Accordingly to the 
formula (4.4) the matrix T can be represented in the form of a product of the matrices Ti t2 ,3 
and Th, each of the same dimension (2 x 2 NM ): 



where 



t = n^nn, (4.7) 

rp a\\...auM T~\ Kza nrn b nrn (a q\ 

J 3,6 U -6 W M _ ll e ' ^-°) 

nm 

CtnmQ>n,rn-\-l ( A Q\ 

- L 2,b 11 ...b NM — u aub 11 ---'J aNM b NM , \^-^J 

nm 

rp ari-..a NM r r TT _Ki 

J l,feu...6jVM — °aii&ii-"°ajVM6]VM ll e ) l^- iu J 

nm 

/ti an...ajVM £ X TT a ha nm { A 

1 h,b 1 i-b N M ~ ailbii-~ a NM bNM 11 e ■ ^■ J "U 



nm 



Here we introduced a new way of indexing the matrix elements in the expression (4.4): 

{auk ■ ■ ■ OjVAffc} = { a ll ■ ■ ■ UNm}, {Cll,fc+1 • ' ' &NM,k+l} = {&11 • ' ' &JVM/, 

and we will continue with these assignments till the end of the paper. 

Further as is known |IB|, if we introduce three sets of 2 NM - dimensional matrices (t^' z ) 



of the form 

t*%> z = 1®1®---®t x > v ' z (iVM -faktors), (4.12) 

where the Pauli matrices r x,y ' z are situated in these products at the nm-th place, the matrices 
7^2,3 and T^, (4.8 — 11) can be rewritten in the form: 

Ti = exp K x T* m T* +liTn , T 2 = exp [K 2 ^ r* m T* m+1 , (4. 13) 
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T 3 = (2 sinh 2K,) NM ' 2 exp K* £ r x nm , (4.14) 



T h = exp 



f/»Ec). (4-15) 

\ nm / 



where the quantities and ifg* are connected by the conditions of the form (2.8), and the 
spin Pauli matrices r^' 2 , (4.12) commute one with each other for different (nm) ^ (n'm'), 
and simultaneously for each given nm these matrices satisfy the standard conditions [JT6] . 



It is easy to see that the matrices T\ i2 ,hi (4-13, 15) commute one with each other, but they 
do not commute with the matrix T 3 , (4.14). Obviously, for (h = 0) we obtain the known 
formulas @ for the matrices Ti t2 ,3, describing the three-dimensional Ising model on a simple 
cubic lattice. To the transition to the 2D Ising model in the interaction constants K% and K 2 
corresponds taking {K\ = 0) or [K 2 = 0) and simultaneously removal of summation over n 
(N = 1) or over m (M = 1) respectively. We obtain this way the standard expressions |||8| 
for the 2D Ising model in external magnetic field, and the operator T\ (4.13) is identically 
equal to the unit operator (Ti = 1) in the first case, and (T 2 = 1) in the second case, 
respectively. 

A bit different situation occurs in the case of transition to the 2D Ising model in the 
interaction constant K 3 . In this case one should take (K 3 = 0, K = 1), i.e. omit summation 
over k. As a result one can arrive at the following formula for the operator T 3 (4.14): 

r 3 ^T 3 (^3 = 0)=n(l + <J. (4.16) 

nm 

Namely this structure of the operator T 3 * enables, finally, effective rebuilding of the magnetic 
operator Th (4.15), as it was shown above on the example of the ID Ising model. In this 
case we can write the expression for the statistical sum for the 2D Ising model in the form: 

Z 2 {h)=Tr{T;T 2 T 1 T h ), (4.17) 

where the matrices Ti^h are defined by the formulas (4.13, 15), and the matrix T 3 * is defined 
by the formula (4.16). The advantage of the representation of the statistical sum (4.17) is, 
in the opinion of the author, in a sense obviously. We will write about this issue additionally 
below. As it will be clear from what is stated further the matrix T 2 TiTh can be conveniently 

1/2 1/2 

written in the form T h T 2 T{F h , where we applied commutativity of the factors, following 
from commutativity of the matrices r 2 m . The statistical sum (4.17) we rewrite in the form: 

Z 2 (h) =Tr{T;T 1 h /2 T 2 T 1 T l h ' 2 ), (4.18) 

1/2 

where the matrix T h is defined by the formula 



Ti Til / 2 

h/2 = T h = ex P 



(V2)E 



mil 



(4.19) 



Below we will use both the expression (4.17) and the representation (4.18), having in mind 
further applications in graph theory (TBI [Tj| . 
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V. TRANSFORMATION OF T-OPERATOR 



A. Introduction of Fermion operators 



Schultz, Mattis and Lieb || showed that the T-matrix in its standard representation 
can be expressed in terms of the second quantization Fermi operators. For this aim they 
applied the known Jordan- Wigner transformations |TP[ which enable expression of the Fermi 
operators (cj„,c m ) for the one-dimensional system by the Pauli operators (r^), ||: 



m— 1 \ / m— 1 

,„ - (>xp | ITT T t T I T m' C m = eX P W 

3=1 / V 3=1 



E 



3 3 ™ 



(5.1) 



As was shown in [12], there is an analogue to the Jordan- Wigner transformations (5.1) 
which generalizes the former to the two-, three-, and <i-dimensional systems. 

For this aim we introduce first the following variables M to the formulas (4.13 — 16): 



r* = -(t z ± ir y ) 

nm 2 nm nm) ' 

which satisfy anticommutation relations for the same site: 

r T + T - \ — i ( T + \2 _ ( T - )2 

I'nm' 'nmj + > \'nmJ \ nml J 

and commutation relations for various sites, 

- 0, (nm) ^ (n'm'). 
Quantities r^ m are often called Pauli operators. The correspondences 



L 'nm' n'm'J 



nm V nm ' nm 



1/2), 



enable to rewrite the expressions for Tx^h and T 3 *, (4.13 — 16) in the form: 



T\ = exp 



1 / A ' nm, 1 ' nm J \ ' n+1 ,m 1 ' n+l,mJ 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



T 2 = exp 



•^■2 ^ X ^nm "I - 'nm)('n,m+l ~^ ^~n,m+l) 



(5.7) 



T ft = exp 



h J2( T nm + T nm) 



(5.8) 



r* = n [1 + (1 - 2T+ m T- m ) 



(5.9) 



As it was mentioned above, the Pauli operators T^ m behave as Fermi operators when 
considered for one site, and as Bose operators when considered for different sites. In order 



11 



to transform to the fermionic representation, i.e. to the Fermi operators in the whole lattice, 
we will introduce an analogue of the Jordan- Wigner transformations (5.1), which will enable 
to express Fermi operators (c^ m , c nm ) by Pauli operators r^ m for the two-dimensional system. 



Namely, there exist in the two-dimensional case two sets of such transformations |L2[ , which 
we represent here in the form: 

(n-l M m-l \ 

*^EE T kl T kl T nl r nl T nmi 

k=l 1=1 1=1 / 



/ n-l M m-l \ 

a nm = exp I in £ ]T r fe/ T fe/ + ^ £ ^nZ T nm> ( 5 - 10 ) 

V k=i 1=1 1=1 / 



and 



(TV m-l n-l \ 

fc=l Z=l fe=l / 



(AT m-l n-l \ 

CT EE T fc/ T M + £ T fcm r fcm W (5.H) 
fc=l i=l fc=l / 

It is easy to show, using formulas (5.3 — 4) that the operators (a nm ,a nm ) and {{3 nm , Pnm) 
are Fermi operators in the whole lattice, i.e. they satisfy anticommutation relations for all 
sites: 

{a{ m ,a nm } + = 1, {alJ 2 = {a nm f = 0; {a{ m , a ] n , m ,} + = ... = 0, (nm) ^ (nW), (5.12) 
and analogously for the /5-operators. There are also inverse transformations: 



r nm = eX P 



n—1 M m— 1 



«L> ( 5 - 13 ) 



fc=l p=l p=l 

etc., which can be easily proved by application of the identities 



exp ill J2 TnnJnm = II ( X ~ 2r nm T nm) = H 
\ nm / nm nm 



' nm' 



from which one can easily derive the equalities 

' vi 'H i ' -fi in vi 'i i i 



nm nm ^nm^nm Hnml^nm 



PlmPnm. (5.14) 



The formulas (5.14) express conditions of local equality of the occupation numbers for a- 
and (3- fermions in one site. Further, as it follows from (5.10 — 11) and (5.13), a- and /3- 
operators are connected by canonical non-linear transformations: 

a\ m = exp(mip nm )ft nm , a nm = exp(i%(p nrn )p nm , 
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Vr. 



N m-1 n—1 M 

E E+E E 

fc=n+l p=l k=l p=m+l 



a lp a k P = [• • -HpPkp, (5-15) 



where the operators (p nm obviously commute with the operators (a£ m , a nm ) and (f3 nm , /3 n m)> 



i.e. 



[V?nm, = ■■• = ■■• = [(finm, Pnm]- = 0. (5.16) 

Commutation relations among a- and (3- operators are more complicated. Namely, as one 
can check by direct calculation that the following commutation relations hold: 

{al im ,(3 nm }+ = {Pi m ,a nm }+ = (-1)^, (5.17) 



a 



a i _rt ot 1 -n (n'<n-l, m'>m+l\ , , 

„m, Pn'm'l- — ■ ■ ■ — [OC nm , P n > m >l - — U, I n' > n + 1 W,' < W, — 1 / ' (,0-loJ 



and 

{aw,/W}+ = . . . = {a] im ,/3j, m ,} + = 0, (5.19) 

in all other cases, where ip nm are defined above (5.15). This way we get rather specific 
structure of commutation relations among a- and (3- operators in the lattice, although this 
structure shows some symmetry. Here is the right place to compare the situation described 
above with the situation we get using the second quantization method. For a system com- 
posed of different particle one introduces the second quantization operators of different kinds 
for different particles. The operators connected to either bosons or fermions satisfy the stan- 
dard commutation relations. As far as the operators for different fermions are concerned, it 
is usually assumed without any proof ||20|| , that within the limits of nonrelativistic theory 



they could be treated formally as commuting or anticommuting. Both assumptions lead 
to the same results when the second quantization method is applied. Nevertheless, in the 
relativistic theory, which allows for transmutations of various particles we should consider 
creation and annihilation operators for different fermions as anticommuting. On the other 
hand, in our case we deal formally with " quasiparticles" of the a- and (3- types underly- 
ing separately the Fermi statistics with commutation relations among particles of different 
types being however dependent on relative position of these "quasiparticles" in the sites of 
lattice. Such a situation, as far as is known to the author, was not present in earlier works 
on application of the second quantization method. 



B. The T 1>2>h - and T 3 * - operators 

Before writing the T- operators (5.6 — 9) in terms of Fermi operators, let as make a few 
remarks. First, the operator T 3 * (5.9) can be expressed in terms of a as well as /3-operators, 
because of (5.14): 



n 



i + (— i) a ™ m 



n 



i+ 



(5.20) 
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where the basic in the Fock representation should be chosen as to be expressed in terms 
of the a- or f3- operators, respectively. Second, the operators Ti 2 ,h we can a l so express in 
terms of either a- or f3- operators. Nevertheless, the operator T 2 we write in terms of the a- 
operators and the operator Ti we write in terms of the (5- operators for reasons which will 
become clear later. 

Now, due to (.10 — 11) we can write the operator T h (5.8) in the form: 



T h = exp 



exp 



h^2i/j nm ((3l m + (3 nm ) 



(5.21) 



where 9 nm is defined as the first factor in (5.13), and ip nm is defined by: 

N m-1 n-1 



exp 



fc=l p=l 



fc=l 



Transformation of the operators is a bit more complicated. Taking into account 
cyclic boundary conditions (4.5), we will write first a sequence of equalities analogous to 
(5.14): 



T N,m T l,m 



-l)^ m /3jv, m /3l,m, T N,m T l,m ~ ( l) Nm Pn,™^!,™' 



and 



T N,m T l,m — ( l) Nm /3jV,m/3l, m , T N,m T l,m ~ ( l) Nm PN,mPl,m, 



(5.22) 



-(-i) Mn <M«n,i ? r+ M r nil = -(-l^a^a^i, 



where 



T n,M T n,l = (- l ) Mna n,M a i,H T n,M T n,l = (-l)""^^^,!, 



M 



N 



9n — ( 1) i M n ^ ' ^nm^nmi 



m=l 



n=l 



(5.23) 



$ m = (-1)"™, N m = -£^ m p nm , (5.24) 



which can be obtained by using the formulas (5.10—13). Therefore we can write the following 
representations for the operators Ti i2 : 



Ti = exp 



r M VN-1 



m=l 



71=1 



(5.25) 



N 



T 2 = exp K 2 



n=l 



M-1 



,m=l 



(5.26) 
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Finally, let us express the operator T 2 in terms of the (3 -operators: 



N 



T 2 =explK 2 J2 

(. n=lLm=l 



M-l 



G = exp 



in al ljn a nrn 



exp 



+ (3n ,m+l 

(5.27) 

where the operators G and x nm are defined by the formulas: 



Xr. 



exp 



N 



n-l 



k=n+l 



k=l 



(5.28) 



and we applied the relations analogous to (5.23), but expressed in terms of the (3- operators. 
The operator S introduced above (5.28) is the operator of the number of particles, which is 
connected with the operators iV and M, (5.24) by relations: 



N 



M 



N 



M 



S=Y,M n =J2N m , G = J] 9n = II 9n 



(5.29) 



n=l 



m=l 



n=l 



m=l 



It is easy to see that the operator G, (5.28) commutes with the operators T\ and T 2 (5.25—27), 
but it does not commute with the operator (5.21), because the following relations are 
satisfied: 



\G, /3 nm } + — 



(5.30) 



Of course, we can also express the operators T% and Th in terms of the a- operators and we 
can write down the formulas if they are necessary. 

It was shown above (4.17) that the statistical sum for the 2D Ising model in external 
magnetic field can be represented by the trace of the operator T, which was expressed here 
by the Fermi second quantization operators. Introducing, as in the one-dimensional case, 
a basic in the occupation numbers representation [2D for the a- and (3- fermions (2 NM 
- dimensional space in the Fock representation), and calculating then appropriate matrix 
elements < l\T\l >, it is easy to see that because of multiplicative character of the operator 
T 3 *, (5.20) all matrix elements, besides the vacuum matrix element < 0|T|0 >, are equal to 
zero. For the vacuum matrix element contribution from the operator T 3 * is equal simply to 
2 NM , and we can write Zii)i) (4.17 — 18) in the following form: 



Z 2 (h) = 2 NM < 0| (1^7)010 >= 2 NM < 0\(T h / 2 T 2 TiT h / 2 )\0 > 
where the vacuum state 10 > is defined in the standard manner 



a, 



|0 >= p nm \o >= 0, 



mm) 



1,2,..., 



N(M), 



(5.31) 



(5.32) 



and operators T 12 ,/i are defined by the formulas (5.21) and (5.25 — 27). Let us stress that 
the vacuum state (5.32) for the a- and (3- fermions can differ among themselves at most by 
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a constant phase factor, which in the given case can always be taken to be equal to unity. 
However, it is no longer true in the case of multiparticle states for the a- and f3- fermions, 
because in this case the essential role begin to play phase factors (— l)* 3 "™, (5.15). As an 
exception serve the one-particle states, for which, as it can be easily found from (5.15), we 
have: 



a 1 ' 10 >= (-lY nm B^ 10 >= /5 f 10 > 



for all(nm). In all other cases the a- and (3- states will differ one from each other by 
their sign which depends on indices (nm) of the corresponding sites. This very fact implies 
main difficulty in the proposed approach to solving the problem under consideration. This 
difficulty can be, however, avoided. 

Let us make two remarks here. It is obvious that the representation (55.31) for the 
statistical sum Z 2 (h) does not depend on the kind of variables (a- or (3- operators) with 
which we introduce the basic in the representation of occupation numbers, because equality 
of local occupation numbers (5.14) holds for the a- and f3- fermions. Further, we expressed 
the operator T 2 in terms of the a- and (3- variables (5.26 — 27), although we will work 
mainly with the expression (5.26). The reason is that in the representation (5.27) for T 2 the 
operators x nm are present (5.28). They are phase factors and it is difficult in practice to 
remove them. Difficulty coming from the presence of these operators is of the same kind, 
which was found by the authors of the paper || who considered the case with external 
magnetic field. Simultaneously, the representation (5.26) for T 2 does not involve the phase 
factors, justifying the choice. Nevertheless, the expression (5.27) for T 2 will be necessary 
in the analysis of the boundary conditions, which play here important role. Analogous 
statement applies to the operator T 1; which we expressed in terms of the a- and (3- variables 
(5.25) and which also does not contain phase factors of the type of Xnm- The essence of our 
approach lies in the structure of the transformations (5.10 — 11), which allows for expression 
of the operators Ti i2 in the form (5.26), which does not contain the phase factors. 

Now, we transform the magnetic operator Th (5.21), or more exactly the ket- vector 
T h \0 >, entering the expression (5.31) for Z 2 {h). Exactly in the sense one should understand 
the equivalence of the two operators T h and T£, acting on the vacuum state |0 >. Below 
we will omit |0 >, as this should not lead to misunderstandings. Analogously, we introduce 
the notation T^'L for the transformed bra- vector < 0|l\/2 and the transformed ket- vector 
10 >, respectively, omitting further bra- and ket- vectors of the vacuum state (< 0|, |0 >). 
Continuing with considerations analogous to these, which gave us the expression (3.7) in the 
one-dimensional case, the operator (5.21) we represent in the form: 









" N 


M M-m 


N 


N-n 


M 


T h = (cosh/i)" M exp 


" / j r-'nm 


exp jet 2 


E 


^ / ^ t @nm@n' 


,m+p 


E 


/ j HnmHn-i-kjm 




nm 




n,n' 


m=l p=l 


n=l 


k=l 


m=l 



(5.33) 



where a = tanh h. Analogously, the operators T h % we write in the form: 



^ /2 = (cosh|)^exp 



t Q-nm 




nm 





h/2 

N M M-m N N-n M 

^ , ^ t ^ , (Xn,m+p(%nm ^ , ^ , ^ , Q^n+fc,m'C^ nm 

n=l m=l p=l ra=l k=l m,m' 



(5.34) 
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T h/2 = ( COsh 



h. 



NM 



cxp 





exp |u 2 


nm 





N M M-m 



N N-n M 



n,n' m=l p=l 



n=l k=l m=l 



(5.35) 



where /i = tanh(/i/2). The operators 7^ and T h T 2 are of rather complicated structure. 
However, they do not contain the phase factors any longer. Substituting the expressions 
(5.33 — 35) to the equalities (5.31), the statistical sum Z 2 (h) can be written in the form: 



Z 2 (/i) = 2 JVM <0|(T 2 T 1 T^)|0>, 



(5.36) 



or 



Z 2 (h) = 2 NM < 0| (77T2T1T; + ^AT^T^B^O >= 2 NM < 0|(E/i + U 2 )\0 >, (5.37) 

where the operators Lq j2 are defined in the obvious way, and the operators and T t * r are 
given by the formulas (5.33 — 35), in which one should omit the factors 



ex P "E^m L exp U^cw , exp U^A 

\ nm / \ nm / \ nm 



t 

nm I 1 



and the operators A and B are of the form: 

A = ®nm, B = Y, Pnm- 



(5.38) 



In derivation of (5.36 — 37) we used the fact that the diagonal matrix elements for the 
product of odd number of Fermi operators is equal to zero, and that the following equalities 
are true 



cxp 



a a nm((3 nm ) 



1 + 



where a is a c-valued function. One should remember also that the operators Ti j2 (5.25 — 
27) contain only bilinear products of the Fermi operators. With this ends the rebuilding 
procedure for the magnetic operator. 



C. The Boundary Conditions 

With the aim of further simplification of the operators Ti j2 we should consider boundary 
conditions for the a- and (5- operators, taking periodic boundary conditions for the Pauli 
operators r^ m (5.2) in both indices (nm) as a starting point. Let us shortly discuss this 
problem here. First, since all terms in T 12 and T£ contain bilinear products of the Fermi 
operators, the following formulas are valid: 

[G, Tx]_ = [G, T 2 ]_ = [G, T*]_ = 0, (5.39) 

which shows that the states with even or odd number of fermions are preserved as well for 
the a- as for the (5- particles. The operator G, entering (5.39), is defined above (5.28). 
Analogously, the following formulas are true: 
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[G, U,]_ = [G, C/ 2 ]_ = 0, 



where the operators are defined above. However, the operators g n and g m (5.24) do not 
commute with the operators T = T 2 TiTj* or U = JJ\ + U 2 (5.37), and this fact implies some 
difficulties. First, let us note that the following equalities hold: 

G=n&=n&»» ^=n^=nv. (5.40) 

n m n m 

where by A^, X 9n and X 9m we denoted eigenvalues of the operators G, g n and g m , equal ±1. 

Let us consider first the case corresponding to the state with even number of fermions 
(Xq = +1). In this case, as one can easily see, we should choose antiperiodic boundary 
conditions for the j3- operators with respect to the second index to (for all n), and antiperiodic 
boundary conditions for the a- operators with respect to the first index n (for all to), i.e. 

Pn,M+l = ~Pn,li Pn,M+l = ~ (n — 1, 2, JV); 
a N+l,m = - a l,m, OlN+l, m = -<*l,m, (m = 1, 2, M). (5.41) 

Then the boundary conditions for the f3- operators with respect to the first index n depend 
on g m , and the boundary conditions for the a- operators with respect to the second index 
m depend on g n . More exactly, it depends on at which step we fix the eigenvalues \g m and 
\ 9n , respectively. The only limitations on the choice of the eigenvalues and corresponding 
boundary conditions give equalities (5.40). The whole freedom in the choice of boundary 
conditions consist of 2 N possible boundary conditions for the a- operators in their second 
index, and 2 M possible boundary conditions for the (3- operators in their first index. De- 
tailed analysis shows that we can without loosing generality choose homogeneous boundary 
conditions, i.e. antiperiodic boundary conditions for the a- operators in their second index 
m, which corresponds to (A^ n = +1) for each n, and antiperiodic boundary conditions for 
the (5- operators in their first index n, which corresponds to (\ gm = +1) for each to, i.e. 

a l,M+l = _a n,l> a n,M+l = -«n,l, Xg n = +1, (n = 1, 2, ....N)] 
@N+l,m = —Pl,mi @N+l,m = —fil,m, ^g m = +1, ijn = 1, 2, M)] 

= II X 9, = (+V N = Ii^ m = (+1) M = +1, (5-42) 

n m 

for each parity of the numbers N and M. Analogously, one can show that in the case of the 
odd states (Xq = —1) the boundary conditions for the a- and f3- operators can be written 
in the form: 

«Ul,m = + a l,m. <XN+l,m = +«l,m, (TO = 1, 2, M)\ 
a l,M+l = + a n,l' «n,M+l = +«n,l, ^g n = ~1, = 1, 2, ....JV); 
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Plr+l,m — +Pl,mi Pn+i,™ — +Pl,m, ^g m — (m — 1,2, M); 

fll,M+l = +$1,1 j Pn,M+l = +Pn,l, (n = 1, 2, N)] 



A d =nv = (-1)^=11 v = (-i) 



-1, 



(5.43) 



for A" and M odd. It is obvious that the constraints on parity of N and M are not important 
here, because we can always choose A" and M in the form (N = 2N' + 1, M — 2M' + 1), 
and then go to infinity with N' and M' independently. 

One can show exactly that the boundary conditions for the a- and (5- operators chosen 
this way are not contradictory, if we take into account simultaneously the conditions of local 
equality of the occupation numbers for the a- and (3- fermions (5.14). As a result we can 
write down the following expressions for the operators T i>2 (5.25 — 26): 



if = exp 



N,M 



Ki E (A 



t 

nm 



n,m=l 



(5.44) 



1 2 



exp 



N,M 



K2 E 



){ a n,m+l + a n,m+l) 



n,m=l 



(5.45) 



where the upper sign (+) corresponds to the states with even numbers of fermions (A<a = +1), 
and the lower sign (— ) corresponds to the states with odd numbers of fermions (A^ = — 1), 
with the appropriate boundary conditions (5.41 — 43). This way, the form of the operators 
Ti j2 for the even and odd states is preserved. What is changing is only the boundary 
conditions. 



VI. THE PARTITION FUNCTION 

In this section we perform all the calculations for the statistical sum written in the form 

(5.36) , and in the end of the section we only give the results for Z 2 (h) written in the form 

(5.37) symmetric in the magnetic operator. This can be done almost automatically, since 
all calculations in the latter case are analogous to the ones given below. 

Collecting all the results derived above, we can write the following expression for the 
statistical sum (5.36): 

Z 2 (h) = (2 cosh h) NM < 0|T|0 >, T = T±T±T* h , (6.1) 

where the operators T£ and T± 2 are defined by the formulas (5.33) and (5.44 — 45). In the 
formula (5.33) one should only omit the factor exp(a J2nm Pl, m ) 5 an d move the constant factor 
(coshh) NM out of the expression < 0|(...)|0 >. Let us remind before the diagonalization of 
the T- operator in (6.1), in which the multiplicative components T^ 2 and are expressed 
by the Fermi operators of the a- and j3- types, that these operators satisfy the mixed 
commutative relations (5.17 — 19). As a result also their Fourier transforms will satisfy in 
general rather complex commutative relations. 
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A. Momentum Representation 



Let us pass now to the momentum representation: 

t _ exp(?7r/4) ^ ung+mp) ct o t = exp(i?r/4) ^ i(nq+mp) t , . 

" m (AM) 1 / 2 j£ qqp ' 1 nm (AM) 1 / 2 j£ lqp ' 1 ; 

where the factor exp(i7r/4) was introduced for convenience. Antiperiodic boundary condi- 
tions give: e tqN = — 1, e tpM = — 1, where 

7T 37T 

qip)= ± A(M)' ± A(M)' ± "--' ^ 
and the periodic boundary conditions give: e iqN = 1, e ipM = 1, where 

271 4:71 

*>=°- ± m)' (6 - 4) 

Substituting (6.2) into (5.33) and (5.44 — 45) we get after straight forward transformations 
the following expressions for the operators T^ 2 1 

If = exp < 2K± J2 [HpVqp + vl-pVq,-p + V-qpV-qp + V-q-pV-q-p) COS q+ 

I 0<<7,p<7T 

(»7-« -P^L + ^-M^J -P + Vq,pV-q ,-p + Vq,-pV-q,p) ™ q] } = H lt{q,p), (6.5) 

0<q,p<7T 

Tf = exp I 2K 2 [(fjpfff + + i-qpi-qp + Z-q ,-pt-q -p) COS P+ 

[ 0<g,p<7r 

0<<7,p<7T 

T* = exp 2 [«(^9)(^-ff,- P ^ + ^!,-p)]+ $ w|= II T±(q,p), (6.7) 
where a(h, q) = tanh 2 /i±±^. 

v ' ^ / sin g 

In the formulas (6.5 — 7) the upper sign (+) corresponds to the case of even states, for 
which one should omit the term in the formula (6.7), and the lower sign (— ) corresponds 
to the case of odd states with respect to the operator of the total number of particles (S). 
The function <&(h) is of arbitrary complicated form and we will not write it down here. We 
only mention that it plays an analogous role to the role played by the function /(g) in the 
one-dimensional case (3.11). We used commutativity of the operators T^ 2 (q,p) and T^(q,p) 
for different (q,p) to write the expressions (6.5 — 7). Namely, 

[Tf{q,plTf{qi,p>)\_ = [T±(q,p),T±(q> ,p')]- = [T^q,p),T^q' ,p')}. = 0, 
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as can be easily verified. The statistical sum (6.1) we rewrite now in the form: 



Zf{h) = (2cosh/i) 7VM <0 



NM 



n Tt&p) 



0<q,p<TT 



II T?(q,p)T±(,p) 

0<<7,p<7T 



|0>, 



(6.8) 



where |0 > - the function of the fermionic vacuum in the space of occupation numbers 
in the momentum representation. This function was denoted in the same way as in the 
"coordinate" representation but this should not lead to any misunderstandings. We used 
also commutativity of the operators T±(q,p) and T^(q,p) for (q,p) ^ (q',p') to write the 
formula (6.8). The operators (£ 9P , £| p ) and (rjqp,Vqp 9 ) satisfy the standard commutation 
relations. On the other hand the commutation relations mixing them are of rather complex 
form, in contrast to the relations in the "coordinate" representation (5.17—19). This, by the 
way, is the cause of the lack of commutativity of the operators T 2 ± (g,p) and T^(q,p)T^(q,p) 
for (q,p) 7^ (q',p')- Now we will maximally simplify the bra-vector < 0|(...) and the ket- 
vector (...)|0 >, which are present in the expression (6.8) for Tf(h), "transferring" the 
corresponding operators through the vacuum state. 

Now, we will consider in some details the case corresponding to the even number of 
fermions (A^ = +1) which means the choice of the antiperiodic boundary conditions (6.3). 
In the end of the paper we will shortly consider the case of the odd states (Xq = — 1) 
to which correspond the periodic boundary conditions (6.4). First, let us note that the 
following equality holds 



£qp£lP X/ ^ 



qp^qp 



(6.9) 



q,p 



q,p 



Further, it is obvious that for fixed (q,p) the quantities T^(q,p) and T^(q,p)Tj^(q,p) are 
represented by the matrices of the size (16 x 16), each of which is considered in its own space 
of states and , respectively. After introduction of the bases, each of which is built of 
16 functions: 



$ = |0 >£, $ 
= |0 > rn * 



q-p; q,p 



q,p 



q-p; q,p 



= £ t £ f <3>n 
S— q,~ pSqp^U) 

V-q,- P Vlp^0, 



(6.10) 
(6.11) 



where $o an d are the functions of Fermi vacuum (which, as was mentioned above, we 
denoted by $ — — |0 >), we obtain after a sequence of transformations the expression 
for the statistical sum Z£(h) in the case of the even states: 



Z+{h) = (2 cosh h) 



NM 



where 



II A l(9,h) 

i 0<g,p<7r 



II Al(p) \ <0|T 2 + T 1 + (/,)|0>, (6.12) 



V 0<<7,p<7T 



f+(h) = exp 
f+ = exp 



B l (q) (rjl q _ p rjl p + r]\ p rf q _ p ) 

0<q,p<ir 
0<q,p<ir 



(6.13) 
(6.14) 
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and the quantities Ai(q, h), ... , B 2 (p) are defined by the formulas: 

Ai (g, h) = cosh 2Ki — sinh 2K\ cos q + a(h) sinh 2K\ sin q 
A 2 (p) = cosh2i^2 — sinh2fT 2 cosp, a(h) = tanh 2 h 



2 1 + cos g 



, . a(/i) [cosh 2Ki + sinh 2Ki cos g] + sinh 2Ki sing 

Ai{q, h) 



sing 
sinh2i^2 snip 



(6.15) 



Analogously, we can after a sequence of transformations, similar to given above, get the 
following expression of the statistical sum understood in the form (5.37) symmetrical with 
respect to the parameter of the external magnetic field h: 

Z 2 + (/>) = (2cosh 2 V A/ ( 1] C 2 M ) ( II Cl(p)) <0\V 2 + (h/2)V+(h/2)\0>, (6.16) 

\0<g,p<7r / \0<q,p<iT 

where 



V+(h/2) = exp 
V+(h/2) = exp 



J2 D 1 (q,h/2)(rji qj _ pV l p + rj 

0<q,p<n 



q,pVq,—p) 



D i (p> h l 2 ) (Up^-q-p + ,-p) 

0<q,p<7r 



(6.17) 
(6.18) 



and the quantities C*i(g, h/2), ... , /^(p, V 2 ) are defined by the formulas: 

Ci(g, h/2) = cosh 2Ki — sinh 2Ki cos g + a(h, q) sinh 2i^i sin g, 

C*2(p, h/2) = cosh2i^ 2 — sinh 27^2 cosp + a(h,p) sinh 2.^2 ship, 
a(h, q) [cosh 2K± + sinh 2K± cos g] + sinh 2Ki sin g 



DM h/2) 
D 2 (p, V2) 



Ci(g,V2) 

a(/z,p)[cosh2fT 2 + sinh 2i^ 2 cos p] + sinh 2fT 2 sin p 
C 2 (p,V2) ' 



(6.19) 



where a(h,x) = tanh 2 (/i/2)(l + cos x) /(sin 2). For the purpose of derivation of the ex- 
pressions (6.16 — 19) we applied the fact for the even states the operator U 2 in (5.37) gives 
vanishing contribution to the statistical sum Z 2 (h). We gave here two representations (6.12) 
and (6.16) for Z 2 (h), because as can be shown they both can be applied in the graph theory 
|15| as we mentioned above. 

In principle, one could now expand the vacuum matrix element in (6.12) or in (6.16) into 
a sum of vacuum matrix elements of the type 



q-p 



' V-q'-p'rfq' 



|0>, 



derive appropriate commutation relations for the £- and r^-operators, and, finally, sum up 
the series. Nevertheless in practise this task seems to be extremely difficult, as believes 
the author. Therefore we will proceed the other way. Namely, we will come back to the 
"coordinate" representation, i.e. to the a- and (3- operators. Then the operators T^ 2 (6.13 — 
14) or V12 (6.17— 18) are expressed as follows: 
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" N N-n M 

v, + = exp EEE <i)0LJin, 

_n=l 1=1 m=l 
' N M M~m 
V 2 + = exp E K k ) a n,m+kOi n 

_n=l m=l k=l 



(6.20) 



where a(l) and b(k) are given by: 

<d = 1 

0<(?<7T 

for the "symmetric" case, and by 



]T 2A(g)sin(Zg), 6(fc) = — £ 2£> 2 (p) sin(fcp), 



(6.21) 



0<p<7T 



C (0 = ^ E 2S 1 ( ? )sin(Z g ), d(k) = — J2 2B 2 (p)sm(kp), 



0<q<n 



M 



(6.22) 



0<p<n 



for the " nonsymmetric" case, where the quantities Bi(q), ... , D 2 (p) are defined above by 
(6.15) and (6.19). Here we used in both cases the same notation V{^ 2 , and further we continue 
with this convention. As can be seen from (6.13 — 14) and (6.17 — 18), the structure of the 
operators T^ 2 in the "coordinate" representation is the same as in the case (6.20). The 
only change concerns the weight factors: a(l) — > c(l) i b(k) d(k). The whole procedure 
used above corresponds to the renormalization of the interaction constants in the former 
expression (5.31) for the statistical sum. We will be exploring this topic more thoroughly 
in the following papers of this series. Moreover, here appears also a delicate problem of the 
boundary conditions, connected with the expressions (6.20). The discussion of this problem 
we also postpone to a future publication. Here we mention only that in the thermodynamic 
limit we can neglect the boundary effects. On the other hand, in the situation at hands it 
is much easier and more convenient to consider the diagram representation for the vacuum 
matrix element < 1 l/f 1- 1 > in the "coordinate" representation than in the "momentum" 
one, which we denote here by S, i.e. 



s =< o|\/ 2 + v; + |o >=< ojcjo > 



(6.23) 



B. The Diagram Representation for S 

Our aim now is to calculate the vacuum matrix element S (6.23) for the sum of products 
of Fermi creation and annihilation operators. The operator G entering (6.23) is a polynomial 
in the variables a(l), b(k), a nm and /3^ m . Since G enters in the (6.23) expectation value form 
< 0|G|0 >, not all terms in the polynomial give a different from zero contribution to the 
matrix element S. Expanding G and substituting the expansion into (6.23), the quantity 
S can be represented in the form of the sum of the vacuum matrix elements Y^v where 
S u is the vacuum matrix element for the z/-th term of the polynomial G. As it follows 
from (6.20), all terms of the polynomial G are products of various pairs b{k)a n ^ m+ ka nm and 
a {l)flnmPn+i,mi which we will call below a- pairs and (3- pairs. Obviously, all the terms in the 
polynomial G with non-equal numbers of the a- and f3- pairs give vanishing contribution. 
Moreover, not all terms in the polynomial G with equal numbers of the a- and j3- pairs 
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will give nonvanishing contribution to S. Namely, the non-zero contribution to S will give 
only these terms with equal numbers of the a- and (3- pairs, in which each annihilation 
operator a nm is paired with the corresponding creation operator /?„/ m / with identical indices 
(n = n',m = m'). In the opposite case this term gives obviously no contribution to S. 

This way we arrive at a diagrammatic representation by noticing that to each vacuum 
matrix element S v we can uniquely assign a set of lines (links), connecting some of the sites 
of the lattice. For example, to the graphs at Fig.l, a) — d) correspond the following matrix 
elements: 

a) a(2)6(3) < 0|a n|TO+ 3a nm /3^ jm+3 /3jj +2im+ 3|0 > ; 
b) a 2 (l)a(2)6 2 (l)6(2) < 0| 

^■Pn+l,m— l@n+2,m— lPnm^n+2,m(^n,m+xPn+X,m+l\^ ^ ) 
c) a 2 (l)a 2 (4)6 2 (l)6 2 (2) < 0|a; n|7n +x«nm a n+l,m+l a n+l,m«n+l,m-2CK n 4.i )m _4a; n +5 )m _2a; n +5 )Tn _4 

vflt at /at /at /at ot q\ /at If) > ■ 

^ r J n+l,m— 4rn+5,m- 4rn+l,m- 2r / n+5,m— 2rnmrn+l,mrn,m+lrn+l,m+l I w ) 
d) a 2 (2)a(4)6(2)6(3)6(5) < 0|a nim+2 a nm a„ + 2, m an+2,m-3«n+4,m+2an+4,m-3 

X /^n+2,m-3/^n+4,m-3/^nm/^n+2,m/^n,m+2/^?i+4,m+2|0 > • (6. 

As one can see from the formulas (6.20) and (6.24), to each horizontal line of the "length" k 
corresponds the factor b(k). Also, to each vertical line of the "length" I corresponds the factor 
a(l). The a(l) and b(k) are defined by the expressions (6.22) for the nonsymmetric case. 
As was shown above a nonzero contribution to S give only these matrix elements S u , which 
do not contain equal numbers of the a- and (3- pairs. Moreover, the necessary condition 
for a non-zero contribution is the annihilation operators a nm pair with the corresponding 
creation operators f3^ m . Geometrically this condition means that from the whole family of 
possible graphs only those for which in each site meet under "right angle" only zero or 2 
lines (links) give a non-zero contribution to S. In other words, the graphs in any site of 
which meet two horizontal or two vertical lines are forbidden. The simplest examples of such 
graphs are shown in Fig.l, b) — e). As a result all graphs giving non-vanishing contribution 
toS should be closed. Moreover, in each site of the graphs selfintersections of lines (links) 
are forbidden, since (a nm ) 2 = (Pl m ) 2 = 0. ^From the point of view of the graph theory to 
the closed graphs described above correspond nonoriented Hamilton cycles (with valency of 
sites 5 = 0, 2) on the simple rectangular lattice [|l8l.|23|. pf. 

This way the vacuum matrix element S (6.23) can be represented in the form 

S = S u = \sum on all closed graphs], (6.25) 

V 

where in the calculations every multiple-connected graph is counted as one (for example, 
the graph in the Fig.l, c)). Every closed graph gives the contribution equal to 

(±l)f[a(l j )b(k j ), (6.26) 

i=i 

where s is the number of the horizontal links, which is equal to the vertical links. Further, 
applying the connection between the a- and (3- operators (5.15 — 19), and the Wick theorem 
[l].[22|], one can show that any vacuum matrix element giving non-zero contribution into 
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the sum S (6.25), can be split into a product of the matrix elements, corresponding to 
the connected parts of the graph (which we will call below for shortness the simple loops 
without selfmtersections in the sites of the lattice). We can check by direct computation, 
using the commutation relations (5.17— 19) for the a- and (3- operators, that for example 
the graphs from the Fig.l, (b — d) contribute with the sign (+). Other graphs can contribute 
with the sign (— ) as well as, for example, the graph in Fig.l, e. Commutation relations for 
the a- and (3- operators (5.17 — 19) are illustrated in an appealing way in the Fig. 2, where 
the distinguished operator a nm (*) for the fixed site {nm) commutes with the (3- operators 
in the sites (n ; m'), signed with the cross (x). For all others sites the a- and (3- operators 
anticommute. As a result the contribution from each particular graph splits into a product of 
contributions from the simple loops. The contribution from a simple loop with s horizontal 
and s vertical links is equal to: 

3=1 

The expression for S (6.25) is now of the form: 

{'} {s},{q} 

where a(lj) and b(kj) are functions of z\ = tanhiTi, z 2 = tanhi^ 2 and y = tanh 2 (/i/2) for 
the symmetric case, and y = tanh h in the case asymmetric with respect to the parameter 
of the external magnetic field (h). A contribution to (6.28) gives besides summation over 
the number of links s, also the summation over all lengths of these links {k} and {I}, for 
fixed s. As can be easily seen, the summation in (6.28) over the lengths of the horizontal 
{k} and vertical {1} links is performed independently. In the graph theory [18] , 2^ the 



function (6.28) is called the generating function, as we mentioned above, introducing for it 
the notation T^(zi, z 2 ,y), where the upper index (h) means being a member of the set of 
Hamilton cycles. The problem was reduced this way to the summation over all Hamilton 
cycles with the varying length of the step (edge) on the rectangular lattice of the type 
described above. 

Now, let us note that the graph representation of Z 2 (h), described above, looks similar 
to the diagrammatic representation for the statistical sum of the 2D Ising model in the 
vanishing magnetic field (h = 0), (see, e.g., the papers p5| - |27| ). In this case, as is known 



25[, the statistical sum can be represented in the form: 



Z(K 1} K 2 ) = (2 cosh Ki cosh K 2 ) NM £ g a3 tanh a K x tanlr 3 K 2 , (6.27) 

where g a $ denotes the number of the closed graphs consisting of (3 horizontal and a vertical 
links. Since these links connect the closest sites of the square lattice, to each link a is 
assigned the factor (weight) tanhi^i, and to each link (3 is assigned the factor tanhi^ 2 . In 
some sites of the graph a simple selfintersection is possible, i.e. in one site of the graph 
meet zero, two, or four lines. This corresponds to the nonoriented Euler cycles of the degree 
5 < 4, |T8"|, |24|1 . In the Fig. 3 is shown one of the simplest graphs contributing to the sum 
(6.29) for Z 2 (Ki, K 2 ). The essential difference of this case in comparison with the case with 
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the field (h), described by us above, lies in the latter property, because in our case in one site 
of the lattice can meet only zero or two lines (horizontal and vertical). This corresponds, 
as was discussed above, to the nonoriented Hamilton cycles on the square lattice [|T8"|, |2"4"f . 
The second difference is that the a- and (3- links in (6.28) can connect not only the nearest 
sites of the lattice. This result in the appearance of dependence of the weight factors a(lj) 
and b(kj) on the distances I and k between the sites of the lattice in the vertical and the 
horizontal direction, respectively. As we mentioned above, the problem of calculation of the 
sum (6.28) can be called in the language of the graph theory |18j the problem of summation 



over the Hamilton cycles (simple cycles) on the rectangular lattice with (N x M) sites with 
varying "length" of the edges in the horizontal and in the vertical directions, respectively. 
Simultaneously, the problem of the sum (6.29) is equivalent to the problem of summation 
over all possible Euler cycles, described above of the type (5 < 4) on the same lattice. As 
is known there is a close correspondence between the Euler and the Hamilton graphs. 



For some types of the Euler graphs one can consider instead the corresponding Hamilton 



graphs. The reversed statement is not true. In the papers [15j is shown one more example 
of the nontrivial connection between the generating functions for the Euler cycles and the 
Hamilton cycles on the simple rectangular lattice. Namely, in the papers Jl5[ was shown 



that the generating function T( h >(zi, z 2 ,y = 0) for the Hamilton cycles described above is 
exactly equal to the generating function r^(z±, z 2 ) for the Euler cycles (5 < 4) for the 2D 
Ising model fl8 |. Therefore, the following equality is true: 



rW(^ 2 ) = n n 



l+zl){l+zl)-2 Zl {l-zl) cos ^—2z 2 {l-zl) cos ^ 



=1 m=l 



(6.28) 



where Zi j2 = tanhi^ 12 . Taking in (6.16) the external magnetic field to be equal to zero 
(h = 0), and using the equality (6.30) we arrive at the classical expression [|J for the free 
energy on one Ising spin in the 2D Ising model. Let us note that the contribution of each 
graph (connected or disconnected), which consists of a set of the Hamilton cycles, can be 
represented in the form of a product of the determinants of the incidence matrices B v : 

(±)I[det\B u \ , 

where u denotes the order of connectedness of the graph under consideration. It is equal 
to the number of the simple loops creating the graph. This way we conclude that for the 
computation of the statistical sum for the 2D Ising model in the external magnetic field 
it is necessary to calculate the generating functions for the Hamilton graphs on the simple 
rectangular lattice of the type described above (see |0|, |28|, [f29]j). 



VII. LIMITING CASES 

A. The Onsager solution 

Let us shortly discuss one of the method of receiving the Onsager solution pE| . Putting 
in Eqs.(6.16) and (6.19) the magnetic field equal to zero (h = 0), the partition function Z 2 
(6.16) takes the form: 
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NM 



Z 2 = 2 NM [(1 - zf)(l - z$]- ±1 * L < 0\T;T*\0 >, (7.1) 

where z 12 = tanh(i^ li2 ), while operators T^ 2 can be written in the "coordinate" represen- 
tation as: 

{N M N-n ~| ( N M M-m \ 

£ Z lPl,mA+l, m \ > T 2 = exp <Y1 S S 4«n,m+A«nfn [ • (7-2) 

n=lra=l J=l J U=lm=l k=l ) 



The Ward-Kac solution [fjQfl , briefly described in plj , contains a topological considera- 
tions. Namely, for a given closed graph (we consider here Euler graphs on a lattice) a factor 
a = exp(z7r/4) is added to a left turn, and a factor a -1 = exp(— in/ 4) to a right turn. Closed 
graphs (i.e. which we want to include) are thus taken into account and forbidden graphs 
are compensated if me follow various paths over these graphs. Full proof of this theorem 



was given by Sherman p2\ . Similar result holds for hamiltonian graphs on a lattice with 
variable length described above which will be shown in simple cases below. However, we 
will follow the methods of [27, [33| in our consideration. 



First of all let us mention that some of hamiltonian loops (e.g. Fig.l, e) contribute with 
minus (— ) sign in formula (6.28) for S. Namely straightforward verification, with the help 
of commutation relations (5.17 — 18) shows that each doubly intersecting link of the one 
shown in Fig.l, e contributes a minus sign to a overall sign of a simple loop (6.27) for all 
admissible diagrams. At the same time each "simple double link" of the type shown in 
Fig.l, f contributes a (+) sign to the overall sign of a simple loop (6.27). All other simple 

loops without "double links" of the shown in Fig.l, b d come with a plus (+) sign in the 

sum (6.28). (Let us note that there is a one to one correspondence between Euler graphs on 
a lattice and hamiltonian graphs with variable step without " double links" , the hamiltonian 
graph may contain one, two or more simple loops. In order to establish this correspondence 
it is necessary to select in the Euler graph all intermediate vertices together with intersecting 
horizontal and vertical links of the Euler graph.) 

It is easy to understand now, that if in expression (6.28) for S all simple loops are taken 
with a (+) sign, all left (and right) turns in a simple loop give a factor a = exp(i7r/4), 
(a -1 = exp(— z7r/4)), than the problem of calculating the sum for S (6.28) is in fact reduced 
to a "random walk" on a lattice with variable step G7L ETT], 133 . In fact, with such a way of 



following simple loops all loops with "double links" cancel (e.g. loops in Fig.l, e and d), as 
it should be. In this way one can follow all the loops with "double links" and verify that 
they cancel each other. Moreover, one can check using various examples, following the same 
reasoning as given m |7|, §§ that if we follow various paths over all hamiltonian loops 



with variable step without "double link" (including relevant weights a and a~ l at each turn) 
than all the allowed diagrams will cancel. One should stress here that such full cancellation 
of forbidden diagrams in every order takes place only in the case of factorizable weights 
(^1,^2) corresponding to step lengths I and k, respectively. 

Returning to our problem and using the results of p7|,[55|, we obtain for S (6.28) the 
following expression: 

00 

S = expl-Y.fr], (7.3) 

where f r - sum over all single loops with length (r = 2s), i.e. consisting of s horizontal and 
s vertical links. Each horizontal line contributes a factor (z^e^/ 2 ), and each vertical line - 



27 



a factor (z[e J</3//2 ), where angle (cp = ±7r/2) corresponds to left or right turn. Introducing 
quantity W r {n,m,v) - sum over all possible paths with number of links equal to (r = 
si + S2) from a given initial point (no, tuq, vq) to a point (n, m, z/), where z/ - auxiliary index, 
corresponding to four directions (1, 2, 3, 4) on a square lattice, we get for f r : 

fr = y Y W r (n ,m ,v> ). (7.4) 
One can easily get the following recursion relations for W r (n, m, v) with (a = exp(?7r/4)): 

N N 

W r+ i(n, m, 1) = + a^ 1 Y, z[W r (n — I, m, 2) + + a Y z[W r (n + I, m, 4); 

1=1 1=1 

M M 

W r+1 (n, m,2)=a^ z k 2 W r {n, m - k, 1) + + a" 1 4 W r(n, m + k, 3) + 0; 

fc=l k=l 

N N 

W r+ i(n, m, 3) = + W r {n-l,m,2) + + a" 1 z[W r (n + l,m, 4); 

M M 

W r+ i(n, m, 4) = a' 1 4W r {n, m-k,l) + + aY 4W r {n, m + k, 3) + 0. (7.5) 

k=l k=l 

The meaning of recursion relations (7.5) is evident. Since the point (n, m, 1) can be reached 
from (n',m,2) and (n",m,4); i.e. from above and from below (direction "1" was chosen 
to be "right"), where n' — n — l,n" — n + I, and / ranges, strictly speaking, from 1 to 
N — 1. However, for large N the summation over / can be extended to N, which was done 
in expression (7.5), because in the thermodynamic limit these boundary conditions do not 
play a role. Hamiltonian structure of simple loops is evident in the structure of recursion 
relations (7.5), which should be compared to the case of Euler graphs |3~l], |3~3|. Writing the 
relations (7.5) in matrix form 

W r+ i (n, m, u) = Y, A(n,m,v\n' ,m' ,v')W r (n ,m' ,1;'), (7.6) 

n',m',i/ 

one can easily see that the following relation holds: 

TrA r = Y W r (n ,m ,u ), (7.7) 

n ,m ,u 

and also 

fr = Y r TrAr = Yr^ X "- ( 7 ' 8 ) 

i 

Taking into account (7.4) and (3.6) we get for 5, (7.3) the following relation: 



i 

where Aj - eigenvalue of the matrix A(n, m, u), (i = 1,2, 4iVM). The matrix A(n, m, v) 
can be easily diagonalized over indices (n, m) with the help of Fourier transformation: 
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N,M 



W r (n,m,v)= J2 e 2 -*r nq+ ^ mp W r (q,p,v). 

q,p=0 



(7.10) 



Inserting (7.10) into (7.5), for fixed (q,p) we get: 
A(q,p,u\q,p,i/) = 



a- 1 EiLi4 £ ~ lq aT,iLi4e lq 

aT,f=i4e- lq or 1 E^i4^ 9 

a' 1 EfcLi 4u~ kp aJ2k =1 zluj k P 



(7.11) 



where a = exp(i7r/4), e = exp(27T2/iV), u; = exp(27ri/M). 

It is evident, that for fixed (q,p) it suffices to calculate the determinant of (4 x 4) matrix: 



11(1 — Xj) = Det(5 pv , - k w ,) = A(q,p), 

3=1 

and after simple calculations for A(q,p), (7.12) we get the following formula: 

(1 + zl)(l + zl) - 2z x {l - z 2 2 ) cos(2nq/N) - 2z 2 (l - z\) cos(2vrp/M) 



A(q,p) 



2zi cos(2Trq/N) + z?)(l - 2z 2 cos(27rp/M) + z\\ 



(7.12) 



(7.13) 



In (7.13) we have neglected the terms proportional to z x and z 2 , since for large N and M 
i z 2 ~ 0, for zi i2 < 1. Finally for asymptotically large (N, M) for S (7.9) with the 



help of (7.13) we get: 



N,M 



N,M 

n 

<j,p=0 



(1 + z?)(l + s|) - 2zi(l - 4) cos(27rg/iV) - 2z 2 {l - z\) cos(2np/M) 
(1 - 2z x cos{27cq/N) + z{){l - 2z 2 cos(27rp/M) + z\) 



1/2 



• (7-14) 



Of course, for asymptotically large (N,M) the expression (7.14) goes to expression (6.30), 
because of following relations 



N 



IJ (1 - 2*i cos ^ + zl) = 1, n (1 - 2 ^2 cos ^ + z 2 2 ) = 1, 

for (N,M — > oo), *x,2 < 1- Finally, using (6.23) and inserting (7.14) into formula (7.1), 
for free energy per Ising spin in the thermodynamic limit we get the well known Onsager 
solution || . The method of finding the Onsager solution, given in this paper, disregarding its 
complications, allows for analytical study of the Ising-Onsager problem in external magnetic 
field in several limiting cases in two and three dimensions. The proposed method of receiving 
the Onsager solution, as well as previously known graphical methods, work only for case 
(a(Z) = z[, b(k) = z 2 , l(k) = 1,2,...). Can be shown that all these methods are not 
applicable if faktors (a(l) and b(k)) have different functional structure. Anyhow, contrary 
to all previously derived methods (graphical et al.) the presented method allows, in such or 
other approximation, for accounting external magnetic field H. 
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B. Low-temperature asymptotic for F 2 D{h) 



The aim of this chapter is to consider the free energy per one Ising spin in the external 
magnetic field for some limit cases. For that reason the parameters (Ki,2, h) are to be 
renormalised in the following way > 0): 

sinh2K* 2 = /3i j2 [sinh 27^2(1 - tanh 2 (/i/2)], 
cosh(2iT 1 * 2 ) = p li2 [cosh2K h2 + tanh 2 (/i/2) sinh 2Ki j2 ] , 

(3 h2 = [1 + 2tanh 2 (V2)sinh2K 1)2 e 2 ^> 2 ]~ 1/2 , tanh 2 h\ 2 = tanh 2 (/i/2) /?1 ' 2 exp j 2Kl ^ (7.15) 

cosh -Fl\ 2 

The above presented formulae are adequate for symmetrical case. For asymmetric case it 
is sufficient to substitute, for instance, K% — > K 2 ,h 2 — > 0. In short, in this case only the 
parameter K\ and the field h are subjects of renormalisation. Formulae (6.21) and (6.22) 
included in (6.20), take the form: 

a(l) = zf + tanh 2 K ^~ Z \ b{k) = zf + tanh 2 h* ]' 4 " (7.16) 

I 1 ~~ Z l ) I 1 ~ Z 2) 

for symmetrical case and 

a(l) = z{ 1 + tanh 2 h* ]~*[ , b{k) = z k 2 , (7.17) 
I 1 ~ z i) 

for asymmetrical one. 

Equations (7.15 — 17) and the way they were derived point on the possibility of obtaining 
series of asymptotics for free energy per one Ising spin for 2D Ising model in the external 
magnetic field (if). In paper Jl5] has been shown that vacuum matrix element S =< 
0\V 2 + Vi~\0 >, appearing in (6.16) for Z 2 (h), for case (a(t) = y, b(k) = z 2 ) is equal to: 

S = T^ h \z 2 , y) = I] [l + 4 + 2z 2V - 2*2(1 - y) cos(p)] 2 , (7.18) 

0<q,p<n 

The above formula may be used to obtain low-temperature asymptotic solution for the 
free energy F 2 u(h) per one Ising spin in the thermodynamic limits. Note that the condition 
[tanh 2 h*/(l-z*) 2 ] -> 1, together with (7.15) is equivalent to (exp(-2^i)(l -tanh 2 h) -> 0). 
For given J\ = const, H = const the above formulated condition is fulfilled for temperature 
area T, when h ~ e^ 1 , e <C 1. For that reason, if for instance ([1 — tanh 2 h* / (1 — z\) 2 } ~ e, 
then a(l) = tanh 2 /r7(l - 4) 2 + ~ ezf. Consequently in this case the result (7.18) may 
be applied. To prove it let us consider Eq. (6.15) for Bi(q,h), expressed by renormalised 
parameters (h*, K*[): 

tanh 2 h* T ^- + 2z 1 *sino 

BAq, h) = L_l (7.19) 

uy ' ; l-2zfcosg + ^ 2 ' v ; 

where z\ = tanhi^, a h* i K\ connected with h i K\ as was shown in (7.15). Moreover, 
due to identity 
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z\ z\ (1 — tanh 2 h) 



1 + zf 1 + 2 Zl tanh 2 h + z\ ' 

introducing a small parameter (1 — tanh/i) ~ e, e C 1, and developing Bi(q, h) into series 
along e (z\ ~ e), we obtain 

„ , , s (tanh 2 /i* + 2zi) sin o 9 

B 1 (q, h) = ± y. ^+ ~ e 2 

1 — cosg 

Substituting the last expression to (6.22) we come to the formula 

a(l) = tanh 2 /i* + 2z 1 *, (7.20) 

describing a(l) with exactness to the second power of e (~ e 2 ), i.e. in this approximation 
a(l) does not depend on /. Finally, substituting in (7.18) y for a(l) expressed by (7.20) we 
receive in the limiting case the following expression for the free energy F 2 D{h): 

- (3F 2D (h) x ln(2 cosh if* coshf^ cosh/i) + ^- T ln[l + z 2 + 2^ 2 (tanh 2 h* + 2z{) - 2z 2 (l - 

2n Jo 

tanh 2 h* - 2z{) cosp]dp, (7.21) 

where h* and K\ depend on h and K\ according to (7.15). Note that the derived approxi- 
mation (7.21) may be also applied to the case of comparably strong magnetic field (H) for 
which (1 — tanh/i) ~ e, £< 1, (T = const). 



C. High-temperature approximation 

In the range of high temperature we impose {J\^jk B T ~ e), e <^ 1, (Ji y2 = 
const, H = const), i.e. z\ t2 = tanhi^i^ ~ e. In this approximation the bra-vector < 0\T 2 , 
expressed in terms of a-operators by (5.45), can be written as: 

N M 

< 0\T 2 ~< 0| exp(2; 2 ^ ^ Pn,m+lPnm), 
n=l m=l 

i.e. expresses in terms of /9-operators, multiplying all phase coefficients ip nm (5.15) by bra- 
vector < 0|. It allows for diagonalisation of the operator T = T 2 T{F^ in (6.1) and calculation 
of the vacuum matrix element < 0|T|0 >. We will not consider the expressions for the free 
energy, as the mentioned above approximation seems to be crude approximation and are 
not of the special interest. 



VIII. CONCLUSIONS 

The case of infinitely small external magnetic field is very interesting {h ~ e, e 
1, T = const). Because in Eqs. (6.15) and (7.15) the magnetic field h appears in tanh 2 h 
function, the computations should be carried out up to the second term (e 2 ) inclusive. 
The presented approach allows for respective calculations, nevertheless they are long and 
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complicated enough to present them in another paper. We should only like to note here 
a case connected with calculations of the free energy for the external magnetic field H 
asymptotic tending to zero, i.e. fulfilling the condition (h — > 0, N,M — > oo). Neglecting 
in Eq. (7.16) for a(l) and b(k) terms proportional to tanh 2 h\ 2 ~ tanh 2 h/2 for a(l) and b(k) 
we obtain the following asymptotic expressions: 

a(l) x zf , b(k) x z*2 , (h — > 0, T = const). 

In this case we can automatically derive the expression for the free energy, substituting in 
(7.14) z 1>2 for 3* 2 : 

-pF 2D (h^0) xln2 + 21n(cosh/i/2) + 

— - / / In [cosh 2K{ cosh 2K 2 — sinh 2^ cos q — sinh 2^ cos pldqdp, 
2n z Jo Jo 

where cosh 2i^ 2 and sinh2i i f* 2 are defined by Eqs. (7.15). This is the leading asymptotic 
term and the latter for (h = 0) given Onsager solution. The procedure is equivalent to 
considering the asymptotically vanishing magnetic field h in the zero-order approximation, 
which in the author's opinion is worth analyzing. 

The above presented approach to the Lenz-Ising-Onsager problem, on the example of ID 
and 2D Ising model in the external magnetic field may be extended on the 3D Ising model in 
the external magnetic field for the purpose of obtaining the low-temperature approximation. 
All calculations are then, in fact, same as the ones leading to Eq. (7.21), apart from details 
connected with dimension of the considered system. The obtained results will be a subject 
of a future paper. 
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FIGURES 



FIG. 1. Examples of some graphs: (a) - a self-avoiding walk; (b)- (f) - Hamilton cycles on 
a rectangular (N x M) lattice with equal numbers of vertices and edges (with varying length of 
steps). 

FIG. 2. "Geometry" of transposition relations for a- and [3- operators: * - a-operator; x - 
/3-operator. 

FIG. 3. The simplest example of a graph (Euler cycles) giving a contribution to the sum over 
states Z(Ki,K 2 ). 
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